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Abstract: We use solar neutrino data to derive stringent bounds on Majorana neutrino 
transition moments (TMs). Such moments, if present, would contribute to the neutrino- 
electron scattering cross section and hence alter the signal observed in Super-Kamiokande. 
Using the latest solar neutrino data, combined with the results of the reactor experiment 
KamLAND, we perform a simultaneous fit of the oscillation parameters and TMs. Further- 
more, we show how the inclusion of data from the reactor experiments Rovno, TEXONO 
and MUNU in our analysis improves significantly the current constraints on TMs. Finally, 
we perform a simulation of the future Borexino experiment and show that it will improve 
the bounds from today's data by order of magnitude. 



1. Introduction 

Present solar ^, |^ |5| and atmospheric neutrino data |p provide the first robust ev- 
idence for neutrino flavour conversion and, consequently, the first solid indication for new 
physics. Neutrino oscillations constitute the most popular explanation for the data (for 
recent analysis see Refs. [0, ^, and are a natural outcome of gauge theories of neutrino 
mass [10, Besides neutrino oscillations, non-zero neutrino masses manifest themselves 
also through non-standard neutrino electromagnetic properties. In the minimal extension 
of the Standard Model (SM) neutrinos get Dirac masses together with magnetic moments 
(MMs) although these MMs are too small to see any effect at current experiments. 
However, if we consider some extensions beyond the SM, neutrinos may have MMs of or- 
der io-^°-io-i^ IXb, relevant for the present sensitivities. In particular, if we assume that 
neutrinos are Majorana particles, then the structure of their electromagnetic properties 
is characterized by a 3 x 3 complex antisymmetric matrix, the so-called Majorana tran- 
sition moment (TM) matrix, that contains MMs as well as electric dipole moments of the 
neutrinos f\l 
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2. Theoretical framework 



In this work (based on Ref. [ij]), we consider that only three Ught neutrinos exist, which is 
well-motivated by recent global fits of all available neutrino oscillation data [^] . We restrict 
our analysis to the LMA-MSW solution of the solar neutrino problem, as indicated by recent 
global fits of solar data 0, |^, |^ and confirmed by the reactor experiment KamLAND |16|. 



On the other hand, we assume that neutrinos are endowed with TMs and have Majorana 



nature, as expected from theory ||ll|, [ij] 



In experiments where the neutrino detection reaction is elastic neutrino-electron scat' 
tering, like in Super-Kamiokande (SK), Borexino and some reactor experiments 19 1 

the electromagnetic cross section is [pOt] 



dfTpm a^vr /I 1 



;F--^)^eff, (2.1) 



dT rnljjL^ \ T E„ 

where the effective MM square is given by |2] 

^xl^ = a^AUa- + a5_AAl'a+ . (2.2) 

The electromagnetic cross section adds to the weak cross section and allows to extract 
information on the TM matrix A. In this cross section, T denotes the kinetic energy of the 
recoil electron and Ey is the energy of the incoming neutrino. According to Eq. ( p.l[ ), one 
notices that the electromagnetic contribution to the total cross section is more important 
at low energies and therefore experiments observing neutrinos with lower energies will be 
more sensitive to TMs. The 3- vectors a_ and a_(_ denote the neutrino amplitudes for negative 
and positive helicities at the detector. The square of the effective MM given in Eq. ( p.2| ) 
is independent of the basis chosen for the neutrino state |21|. In what follows we consider 



both the flavour basis and the mass eigenstate basis. We will use the convention that a^^ 
and A = (A^/j) denote the quantities in the flavour basis, whereas in the mass basis we will 
use dzf and A = (Ajfc). The two bases are connected via the neutrino mixing matrix U : 



[/^a_ , a+ = U^a+ , A = \U . (2.3) 



Taking into account the antisymmetry of the transition moment matrix A for Majorana 
neutrinos, it is useful to define vectors A = (A^) and A = (Aj) in the flavour and mass 
basis, respectively, by 

= £al3-y-^-f and Xjk = £jkl^l ■ (2.4) 

Thus, in the flavour basis we have Ag^ = A,-, A^t = Ag and Are = ^^i- Note also that 

|A|2 = ^Tr (^A^a) |A| = |A|. (2.5) 

This means that, if we are able to find a bound on |A|, we have not only constrained the 
TMs in the flavour basis but also in the mass basis. 

Now we will present the form that the effective MM square fj.'^g in Eq. (^]^) takes in the 
cases of solar and reactor neutrino experiments. The detailed derivation of the following 
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expressions can be found in Ref. Q. In the context of the LMA-MSW solution of the 
solar neutrino problem, we obtain the effective MM square 

I^Ima = |AP - IA2P + P^r (IA2P - lAiP) . (2.6) 

Here Aj are the components of the TM matrix in the mass basis, and P^i corresponds to 
the probability that an electron neutrino produced in the core of the sun arrives at the 
detector as the mass eigenstate ui in a 2-neutrino scheme. 

In the case of reactor neutrinos, the //^g relevant in reactor experiments is given as 

^i=|A^P + |A.|2. (2.7) 

From this relation it is clear that reactor data on its own cannot constrain all TMs contained 
in A, since Ag does not enter in Eq. (|2.7D . In order to combine reactor and solar data it is 
useful to rewrite /U^ in terms of the mass basis quantities 

^1= |A|2-c2|Aip-s2|A2p-2sc|Ai||A2lcos(5, (2.8) 

where c = cos 9 and s = sinO, 9 being the solar mixing angle. Note that the relative phase 
6 = arg(A]^A2) between Ai and A2 appears in addition to |A|, |Ai| and IA2I. 

3. Statistical Method 

In the following we will use data from solar and reactor neutrino experiments to constrain 
neutrino TMs. The x^-function obtained from the data depends on the solar oscillation 
parameters tan^ 9 and Am^ as well as on the elements of the TM matrix A. Regarding the 
dependence on the oscillation parameters we will take two different attitudes. One is to 
assume that tan^ 9 and Am^ will be determined with good accuracy at the KamLAND 
experiment, and hence, we will consider the fixed points in the tan^ 9 — Arn^ plane 

(method I). In the second approach we will derive a bound on the TMs by minimizing 
with respect to tan^ 9 and Am^ (method II). This second procedure takes into account the 
present uncertainty of our knowledge of the oscillation parameters. 

Let us describe in detail how we calculate a bound on the TMs. Our aim is to constrain 
|A|, therefore it is convenient to consider the as a function of the independent parameters 
|A|, |Ai|, [A2I and 6. The x'^-functions which we are using for the individual data sets (solar 
rates, SK recoil electron spectrum, reactor data) will be described in detail in the following 
sections. When performing the fit to the data, we find that in general the minimum of the 
occurs close or outside the physical boundary of the parameters |A|, |Ai|, IA2I. To take 
this into account we apply Bayesian techniques to calculate an upper bound on |A|. We 
minimize the with respect to |Ai|, IA2I and 6 for each value of |A|, taking into account 
the allowed region < |Ai|2 + IA2P < |A|2: 

x'(|A|)=Min[x'(|A|,IAi|,IA2U)] . (3.1) 

In method I we do this for fixed values of the oscillation parameters, scanning over the 
LMA-MSW region, whereas in method II we minimize also with respect to tan^ 9 and 
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Figure 1: Contours of constant Pel- The shaded region is relevant for Borcxino, whereas the region 
below the dash-dotted hne is relevant for SK. 



Am in Eq. (3.1). Then the x is transformed into a likelihood function 



via 



£ oc exp 



(3.2) 



Now we use Bayes' theorem and a flat prior distribution p(|A|) in the physically allowed 
region, p(|A|) = 0(|A|), to obtain a probability distribution for |A|: 

£(|A|)G(|A|)d|A| 



/(lADdlAI- ^c.^^,^,,^^,^,, . 
An upper bound on |A| at a C.L. a is given by the equation 



/(|A|)(i|A| 



a . 



(3.3) 



(3.4) 



Analyzing the minimization of Eq. ( |3.lD in more detail (see Ref. ||l^), we realize that 
the bound on |A| is strongest if Pel = 0.5 because in this situation ^lma (Eil) 
is maximal. In Fig. |^ we show contours of constant Pgi in the tan^ ^ — /^m?/Ey plane. 
For definiteness, the probabilities in the figure are obtained by performing the averaging 
over the production distribution inside the sun for the ^Be fiux most relevant for Borexino. 
However, the probabilities for the '^B flux relevant for SK are very similar. 

The shaded region in Fig. |l] is the region relevant for Borexino, assuming a mass- 
squared difference in the range 10~^ eV^ < Ann? < 10~^ eV^, whereas for SK the region 
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below the dash-dotted hne is most important, due to the higher energy of the *B neutrinos. 
We can read off from the figure that in most part of the SK region Pgi is very smah, which 
means that the sensitivity of SK for |A| is hmited. In contrast, we expect a much better 
sensitivity of Borexino because, in a large part of the relevant parameter space in this case. 
Pel is close to the optimal value of 0.5. 



4. Analysis of solar and reactor neutrino data 

In this section we describe briefly our analysis of solar neutrino data and the data which we 
are using from reactor neutrino experiments. If we neglect the electromagnetic contribution 
to the small elastic scattering signal in the Sudbury Neutrino Observatory (SNO), the only 
solar neutrino experiment whose signal will be affected by neutrino TMs is SK. However, 
the uncertainty in the flux leads to correlations between SK and all other experiments. 
Therefore, also other experiments give some information on TMs and it is important to 



include them in the analysis |22]. 



We divide the solar neutrino data into the total rates observed in all experiments and 
the SK recoil electron spectrum (with free overall normalization) . Subsequently we consider 
the reactor data in combination with the global sample of solar data. 

4.1 Solar rates 

We include in our fit the event rates measured in the chlorine experiment Homestake Q, 
the gallium experiments Sage |2|, Gallex and GNO [^, the total event rate of SK Q based 
on the 1496 days data sample and the result of SNO HI for the charged current and neutral 
current solar neutrino rates. For the analysis we use the x^-function 

xLcs = J2iRj - D,)ivjtTHRk - Dk) . (4.1) 

Here the indices j, k run over the 6 solar neutrino rates, Dj are the experimental rates and 
Rj are the theoretical predictions, which are calculated as described in Ref. [P, except for 
the SK experiment, whose rate includes an extra contribution from the electromagnetic 
scattering, as shown in Ref. |14]. The covariance matrix VJ^^^ takes into account the 



experimental errors and theoretical uncertainties from detection cross sections and SSM 
predictions of the neutrino fluxes (for details see Refs. |^, |23[|). 

4.2 The Super-Kamiokande recoil electron spectrum 

In this section we consider the shape of the SK recoil electron spectrum. The electro- 
magnetic contribution to the neutrino-electron scattering cross section leads to a different 
spectrum of the scattered electrons than expected from the SM weak interaction. There- 



fore, the SK spectral data provide a useful tool to constrain neutrino TMs [24|. 

We perform a fit to the latest 1496 days SK data presented in Ref. Q as 44 data points 
Di, using the following x^-function 

44 

X^pect = E - Di){V^'"'y\aR, - D,) . (4.2) 
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The theoretical prediction Ri is calculated as explained in Ref. |14]. The covariance matrix 



-^^spect pQj^i^g^jj^g statistical and systematic experimental errors, extracted from Ref. The 
is minimized with respect to the normalization factor a in order to isolate only the 
shape of the spectrum. 

4.3 Reactor data 

Data from neutrino electron scattering at nuclear reactor experiments can be used to 
constrain the combination of TMs given in Eq. ( |2.8| ). In fact, these experiments are the most 
sensitive probe for laboratory searches of neutrino magnetic moment because of the high 
flux of low-energy antineutrinos emitted by a nuclear reactor and the good experimental 
control through the on-off reactor comparison. Here we use data from the Rovno nuclear 



power plant |17], the TEXONO experiment at the Kuo-Sheng Power Plant in Taiwan |18| 



and from the MUNU experiment at the Bugey reactor [19|. To include this information in 



our analysis we make the following ansatz for the x^-function 

The sum is over the three experiments. -/V*]-,s is the observed number of events with the 
one standard deviation error o"*, -/V^eak number of events expected in the case of no 

neutrino TMs (only the standard weak interaction) and N^^i^ii^ is the number of events 
due to the electromagnetic scattering of neutrinos with an effective MM /^r. 



5. Bounds on |A| from solar and reactor data 

First we discuss the bounds on |A| from solar data alone. Fixing the oscillation parameters 
at the current best fit point tan^ Q = 0.43, Am? = 6.9 x 10~^ eV^ from the global analysis 
of solar neutrino data including the KamLAND results 0, we obtain with the Bayesian 
methods described in Section ^ the 90% C.L. bound 

|A| < 3.7 X IQ^^^jjLB (best fit point, solar data). (5.1) 

However, such a bound substantially depends on the values of the neutrino oscillation 
parameters. In Fig. |2| we show contours of the 90% C.L. bound on | A| in the tan^ 9 — Am'^ 
plane. We find that the bound gets weaker for small values of Am^, whereas in the upper 
left part of the LMA-MSW region a bound of the order 2 x 10^^'^ /i^ is obtained. 

By combining solar and reactor data we obtain considerably stronger bounds. At the 
best fit point we get at 90% C.L. 

|A| < 1.7 X 10^^*^^B (best fit point, solar + reactor data). (5-2) 

In Fig. m we show the contours of the bound in the tan^ 6 — /S.m? plane for the com- 
bination of solar and reactor data. By comparing Fig. || with Fig. |3| we find that reactor 
data drastically improve the bound for low Am^ values. 
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Figure 2: Contours of the 90% C.L. bound on |A| in units of 10^^°/iB from solar data. The gray 
(hght) shaded region is the Scr LMA-MSW region obtained in the global analysis of solar neutrino 
data (best fit point marked with a triangle), whereas the green (dark) one corresponds to the Scr 
region obtained after the inclusion of the KamLAND results in the analysis (best fit point marked 
with a star), both from Ref. |7|. The dashed line shows the bound obtained at the global best fit 
point. 



In the upper parts of the LMA-MSW region, the solar data alone give already a strong 
bound on |A|. This is due to the fact that there the probability Pel relevant in SK is close 
to the optimal value of 0.5 (see Fig. |l|). In this parameter region the combination with 
reactor data does not improve the bound significantly. In contrast, for low Am^ values the 
probability Pe\ is very small at the neutrino energies relevant for SK, as seen in Fig. |l|. 
Therefore, solar data give a very weak constraint and only the combination with reactor 
data improves the bound. 

Up to now we have calculated bounds on neutrino TMs for fixed values of the oscillation 
parameters iw? Q and Am^ (method I described in Section ^). These results will be 
especially useful after KamLAND will have determined the oscillation parameters with 
good accuracy. In the following we change our strategy and minimize the for each value 
of |A| with respect to tan^ Q and Am^ (method II). This will lead to a bound on |A| taking 
into account the current knowledge concerning the oscillation parameters. To this end 
we make use of the global solar neutrino data and the latest KamLAND results [^], as 
described in Ref. |25] in order to obtain the correct x'^ behavior as a function of tan^ Q and 
Am^. Performing this analysis we obtain the following bounds at 90% C.L.: 



4.0 X 10 [iB (unconstrained, solar + KamLAND data) 

1.8 X 10^^'^/iB (unconstrained, solar + KamLAND + reactor data). 
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Figure 3: Same as Fig. ||, but iiere tiie contours refer to tlie 90% C.L. bound on |A| in units of 
IO^^'^Hb from combined solar and reactor data. 



Together with Figs. || and ^, the bounds given in Eq. ( |5.3| ) show that current solar 
neutrino data can be used to constrain all elements of the Majorana neutrino TM matrix. A 
combination with data from reactor experiments significantly improves the bound on [A|. 



6. Comparison with other works 

In this section we compare the quality of our analysis with respect to two previous works 
doing similar calculations. In particular, we have chosen the first analysis trying to con- 



strain neutrino magnetic moments with solar data, done by Beacom and Vogel [24|, and a 
more recent article presented by Joshipura and Mohanty [p2|. 



In Ref. [24 1, the SK recoil electron spectrum (825 days) is used to constrain the only 
component of the Majorana TM matrix present in a 2-neutrino scheme, A3. They obtained 
IA3I < l.5xW~^^ fiB at 90% C.L. Using the same assumptions, and just including the SK 
spectrum shape (1496 days) in the analysis, we get at 90% C.L. IA3I < lAxlO^^^ i_ib, 
which shows the equivalence of the two methods. Taking advantage of the global analysis 
presented in this work, using all the available information, a better bound is obtained: IA3I 
< 7.8xlO-";UB at 90% C.L. 

In the second work the authors assume the general case of three neutrinos and 
calculate limits on the three elements different from zero of the Majorana TM matrix. To 
do this, they consider just one component non-zero at a time, and obtain a constraint for 
each TM separately in the mass basis. As data sample they use only the solar neutrino 
rates. To compare with their results, we consider just one TM at a time, but we include all 
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Bounds at 90% C.L. 




Beacom & Vogel p4| 


Joshipura & Mohanty 


present analysis 


|Ai| < 
IA2I < 
IA3I < 


l.SxlO-^Vs 


4.3xlO-iVij 

20.5x10-iVb 
3.8xlO-^°/XB 


1.2xlO"iVi? 
I.OxIO-iVb 
7.8xlO~iVi? 



Table 1: Comparison of tlie results from our work and two previous analysis. 



the data sample in our analysis. Results are summarized in Tab. [T|, where we can see the 
improvement of the present analysis with respect the two works considered. 

Finally, we want to remark that our results are the more general presented up to now, 
because they are calculated in a 3-neutrino framework, in a basis-independent way and 
our limits apply to all components of the TM matrix simultaneously. In addition, we have 
used all the available experimental data, and our bounds are stronger than the obtained 
in previous calculations. 



7. Sensitivity on |A| of the Borexino experiment 



Here we investigate the sensitivity of the Borexino experiment |26] to neutrino TMs. This 
experiment is mainly sensitive to the solar ''Be neutrino flux, which will be measured 
by elastic neutrino-electron scattering. Therefore, Borexino is similar to SK, the main 
difference is the mono-energetic line of the ^Be neutrinos, with an energy of 0.862 MeV, 
which is roughly one order of magnitude smaller than the energies of the ^B neutrino flux 
relevant in SK. 

To estimate the sensitivity of Borexino we consider the following x^-function 

A^bins 

xLexino = E (^^' - A)(l^.^^^)-^(iVf - D,) . (7.1) 

Here A'^*'^ is the theoretical prediction for the number of events in the electron recoil energy 
bin i, given by the sum of events from weak scattering, electromagnetic scattering and 
background 

^th ^ ^weak ^ ^em ^^^^^^ ^ ^bg _ (7^2) 

Di is the (hypothetical) observed number of events that, in order to estimate the sensitivity 
of Borexino for neutrino TMs, we assume generated by neutrinos without TMs 

D, = iVf'"^'^ + N^^ . (7.3) 

Hence, we obtain 

i,j=l 

The minimum of this occurs for /Ulma = and is always zero. With this we calculate a 
bound on |A| at a given C.L. That bound corresponds to the maximum allowed value of |A| 



-9- 



International Workshop on Astroparticle and High Energy Physics 



M. A. Tortola 




0.2 0.4 0.6 0.8 

tan^e 

Figure 4: Contours of the 90% C.L. bound on |A| after 3 years of Borexino data-taking in units 
of 10^^°/is. The gray (Ught) shaded region is the allowed LMA-MSW region at 3cr obtained in the 
global solar analysis (best fit point marked with a triangle), and the green (dark) one corresponds 
to the 3ct region obtained including KamLAND in our analysis (best fit point marked with a star) , 
from Ref. Q. The dashed line corresponds to Pei = 0.5 for ^Be neutrinos, and shows the strongest 
attainable limit. 

which cannot be distinguished from |A| = 0, and is therefore a measure for the obtainable 
sensitivity to |A[ at Borexino. The definition of the covariance matrix in Eq. and 
further details about the simulation of Borexino can be found in Ref. jl^ . 

Using the statistical method described previously we obtain for the current best fit 
point tan^ 6 = 0.43, Am^ = 6.9 x 10"^ eV^ the upper bound (sensitivity) 

|A| < 0.28 X lO^^Vs at 90% C.L. (7.5) 

after three years of Borexino data taking. This bound is about one order of magnitude 
stronger than the bound from existing data, which is partly due to the lower energies 
observed in Borexino, in comparison with SK. We have checked that a combined analysis 
of Borexino with existing data (solar and reactor data) does not improve the bound of 
Eq. ( |7.5| ). In Fig. ^we show contours of the 90% C.L. bound in the tan^ 6 — Am'^ plane. The 
strongest attainable limit is roughly 0.24 x 10^^^ fj,B and, in agreement with the discussion 
we gave m Sec. i it is obtained when Pgi = 0.5, as shown by the dashed line in Fig. 
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8. Conclusions 

In this paper we have presented stringent bounds on electromagnetic Majorana transition 
moments (TMs). Such TMs can contribute to the elastic neutrino-electron scattering signal in 
solar neutrino experiments like Super-Kamiokande or the upcoming Borexino experiment, 
as well as in reactor experiments that detect neutrinos through the neutrino-electron elastic 
scattering process. Using most recent global solar neutrino data we have performed a fit in 
terms of the oscillation parameters and the elements of the complex TM matrix A of three 
active Majorana neutrinos. Taking into account the antisymmetry of the TM matrix, we 
have shown that solar neutrino data allow to constrain all Majorana TMs simultaneously in 
a basis independent way, through the intrinsic neutrino property |A|. 

A fit to the global solar neutrino data leads to the bound |A| < 4.0 x lO^^^fXB at 90% 
C.L. We have also considered the role of reactor neutrino data on neutrino TMs, shown 
to be complementary to solar neutrino data. A combined fit of reactor and solar data 
leads to significantly improved bounds: at 90% C.L. we get |A| < 1.8 x W^^^fiB- In 
the near future a precise determination of the oscillation parameters might be possible at 
the KamLAND experiment, which motivated us to scan the tan^ 6 — Am^ plane and to 
calculate the corresponding bound on |A| in each point. These results are shown in Figs. § 
and ^ for solar data and solar + reactor data, respectively. Comparing our results with 
previous calculations done by other groups we have shown that we have obtained the most 
general and stringent bounds presented up to this moment. 

We have also investigated the potential of the upcoming neutrino-electron scattering 
solar experiment Borexino to constrain neutrino TMs. Performing a detailed simulation 
of the experiment we find that it will improve the bound on |A| by about one order of 
magnitude with respect to present bounds. 
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